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Abstract 

The purpose of the "bootstrap program" is to construct integrable quantum field 
theories in 1+1 dimensions in terms of their Wightman functions explicitly. As an 
input the integrability and general assumptions of local quantum field theories are 
used. The object is to be achieved in tree steps: 1) The S-matrix is obtained using 
a qualitative knowledge of the particle spectrum and the Yang-Baxter equations. 
2) Matrix elements of local operators are calculated by means of the "form factor 
program" using the S-matrix as an input. 3) The Wightman functions are calculated 
by taking sums over intermediate states. The first step has been performed for a 
large number of models and also the second one for several models. The third step 
is unsolved up to now. Here the program is illustrated in terms of the sine-Gordon 
model alias the massive Thirring model. Exploiting the "off-shell" Bethe Ansatz 
we propose general formulae for form factors. For example the n-particle matrix 
element for all higher currents are given and in particular all eigenvalues of the 
higher conserved charges are calculated. Furthermore quantum operator equations 
are obtained in terms of their matrix elements, in particular the quantum sine- 
Gordon field equation. Exact expressions for the finite wave function and mass 
renormalization constants are calculated. 
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1 Introduction 



More than fifty years ago, Heisenberg pointed out tlie importance of studying analytic 
continuations of scattering amplitudes into the complex momentum plane. The first 
concrete investigations in this direction were carried out by Jost and Bargmann 0], 
initially for non-relativistic scattering processes. The original ideas turned out to be very 
fruitful and lead to interesting results on-shell, i.e. for the S-matrix as well as off- 
shell, that is for the two-particle form factors P]. Once one restricts ones attention to 
1-1-1 dimensional integrable theories, the n-particle scattering matrix factories into two 
particle S-matrices and the approach, now usually referred to as the 'bootstrap program', 
reveals its full strength. 

The 'bootstrap program 0. for integrable quantum field theories in 1-1-1- dimensions' 
does not start with any classical Lagrangian, but this program classifies integrable quan- 
tum field theoretic models and in addition it provides their explicit exact solutions in 
term of all Wightman functions. We have contact with the classical models only, when at 
the end we compare our exact results with Feynman graph expansions which are based 
on the Lagrangians. 

One of the authors (M.K.) et al. [0. formulated the on-shell program i.e. the ex- 
act determination of the scattering matrix using the Yang-Baxter equations. Off-shell 
quantities, namely form factors were first investigated by Vergeles and Gryanik |Q in the 
sinh-Gordon model and by Weisz in the sine-Gordon model. The concept of general- 
ized form factors was introduced by one of the authors (M.K.) et al. ITDI. In this article 
consistency equations were formulated which are expected to be satisfied by these objects. 
Thereafter this approach was developed further and studied in the context of several ex- 
plicit models by Smirnov ||lT| who proposed the form factor equations (i) — [v) (see below) 
as extensions of similar formulae in the original article |jlO|. The formulae were proven 
by the authors et al. |]12|. There is a large number of more recent papers |T^-|^ on 
form factors. Also there is a nice application |]50| , pl| of form factors in condensed matter 
physics. The one dimensional Mott insulators can be described in terms of the quantum 
sine- Gordon model. 

Finally the Wightman functions are obtained by taking sums and integrals over inter- 
mediate states. The explicit evaluation[| of all these integrals and sums remains an open 
challenge for almost all theories, except the Ising model [^, |53|. A progress towards a 
solution of this problem has recently been achieved by Korepin et al. ||54|| . 

An entirely different method, the Bethe Ansatz was initially formulated in order 
to solve the eigenvalue problem for certain integrable Hamiltonians. The approach has 
found applications in the context of numerous models and has led to a detailed study 
of various mass spectra and S-matrices. The original techniques have been refined into 
several directions, of which in particular the so-called "off-shell" Bethe Ansatz, which 

^Of course one may also adopt a very practical point of view and resort to the well-known fact that the 
first terms of the series expansion of correlation functions in terms of form factors decrease very rapidly. 
Consequently correlations functions may be approximated very often quite well by simply including the 
low particle number form factor into the expansion. 
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was originally formulated by one of the authors (H.B.) ||56| , |57[] , will be exploited for our 
purposes. This version of the Bethe Ansatz paves the way to extend the approach to the 
off-shell physics and opens up the intriguing possibility to merge the two methods, that is 
the form factor approach and the Bethe Ansatz. The basis for this opportunity lies in the 
observation, that the "off-shell" Bethe Ansatz captures the vectorial structure of Watson's 
equations. These are matrix difference equations giving rise to a matrix Riemann-Hilbert 
problem which is solved by an "off-shell" Bethe Ansatz. 



2 The "bootstrap program" 

The 'bootstrap program for integrable quantum field theories in 1+1-dimensions' classifies 
integrable quantum field theoretic models and in addition it provides their explicit exact 
solutions in term of all Wightman functions. The results are obtained in three steps: 

1. The S-matrix is calculated by means of general properties such as unitarity and 
crossing, the Yang-Baxter equations (which are a consequence of integrability) and 
the additional assumption of 'maximal analyticity'. This means that the two-particle 
S-matrix is an analytic function in the physical plane (of the Mandelstam variable 
(pi + ^2)^) and possesses only those poles there which are of physical origin. The 
only input which depends on the model is the assumption of a particle spectrum. 
Usually it belongs to representations of a symmetry]^. 

2. Generalized form factors which are matrix elements of local operators 

°"*(p:,,...,p;io(x)ipi,...,p.r 

are calculated by means of the S-matrix. More precisely, the equations [i) — (v) 
given below are solved. These equations follow from LSZ- assumptions and again 
the additional assumption of 'maximal analyticity'. 

3. The Wightman functions are obtained by inserting a complete set of intermediate 
states. In particular the two point function for a hermitian operator 0{x) reads 



(o|OMO,o),o)^gi/.../^^ 

X {0\O{0)\p,,... ,Pn) 

Up to now a proof that these sums converge does not exists. 
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g -ix Pi 



^Typically there is a correspondence of fundamental representations with multiplets of particles. 



3 



3 Integrability 



Integrability in (quantum) field theories means that there exits oo-many conservation laws 

a^J^(t,x)=0 (L = ±l,±3,...). 

A consequence of such conservation laws in 1+1 dimensions is that the n-particle S-matrix 
is a product of 2-particle S-matrices 

i<j 

If there exist backward scattering the 2-particle S-matrices will not commute and one has 
to specify the order. In particular for the 3-particle S-matrix there are two possibilities 

S^^^ = S12S13S23 = S23S13S12 




N3 1/1 




which yield the "Yang-Baxter Equation" . 

Examples of integrable models in 1+1 - dimensions are the sine-Gordon model defined 
by the classical field equation 



a 



(p{t, x) - (p"{t, x) + — sin P(p{t, x) 







and the massive Thirring model defined by the classical Lagrangian 



Coleman [£8| proved that both models are equivalent on the quantum level. 

Further integrable models are: ZAr-Ising models, nonlinear cx-models, Gross-Neveu 
models, Toda models etc. In the following most of the formulae and explicit solutions 
are given for the sine-Gordon alias massive Thirring model although often corresponding 
results exist also for other models. 



4 The S-matrix 

For the Sine-Gordon alias massive Thiring model the particle spectrum consists of: soliton, 
anti-soliton and breathers (as soliton anti-soliton bound states). Since backward scattering 
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can only appear for particles with the same mass, the two-particle S-matrix is of the form 



u 



t r 
r t 



u 



's6 



\ 



I 



where the rapidity difference Q — |^i — ^2! is defined by 
Pi = mi (cosh Qi . sinh Qj). 

We start with the soliton (anti-soliton) S-matrix: 

C/^'o' . QSS „, QSS J. CSS 

'^afi ~ /\ ■ '->ss — "'-i ^ss—'^-i '->ss — 

a (5 

s — soliton, s — anti-soliton. As input conditions we have: 



1. Unitarity S{-e)S{e) = 1 

u(-e)u{e) = 1 

t{-e)t{e) + r{-e)r{e) = i 

t{-e)r{e) + r{-e)t{e) = o 



2. 'Crossing' 

u(i7r -e)^ t{e) , r{m -9)^ r{e) 

3. Yang-Baxter 

r(^12)M(^13)r(^23) +t(^12)r(^13)t(^23) = «(^i2)r(^i3)M(^23) 



4. 'Maximal analyticity': 

S{9) is meromorphic and all poles in the 'physical strip' < Im^ < tt have a 
physical interpretation; in particular all bound states correspond to simple poles. 
An S-matrix satisfying this condition is also called 'minimal'. For couplings (? < 
(in the language of the massive Thirring model) there are no soliton anti-soliton 
bound states. Therefore in this region of the coupling constant S{9) is holomorphic 
in the 'physical strip' < Im ^ < tt. 



5 



The S-matrix bootstrap using the Yang-Baxter relations was proposed by Karowski, 
Thun, Truong and Weisz 0. It was shown in this article that the 'minimal' general 
solution of these equations is 



, I dt smhf(l - z/) . ^ e 

uW,v) = exp / . , , —r smht— 

^ ^ ' t sinh 4 cosh ^ in 



This S-matrix was first obtained by Zamolodchikov ^9) from the extrapolation of semi- 



classical expressions. It has been checked in perturbation theory. The free parameter v 
is related to the coupling constants by 

1 Stt 2g 
- = 1 = 1 + — . 

The second equation is due to Coleman and the first one is obtained by analyzing the 
pole structure of the amplitude u{9,v). The assumption of 'maximal analyticity' and 
comparison with the known semi-classical bound state spectrum implies the identification 
of the parameter v (see below). 

Bound states 

Let 7 be a bound state of the particles a and {3 with mass 



= -^J m?^ + rn?p + 2mafrip cos a , (0 < a < tt) 
then the 2 particle S-matrix has a pole such that 

^Res5f;">) = rl, 

ti=ia 

P' a' P' «' 

.ResX - X 

a /3 a (3 

where a is called the fusion angle and F^^ is the 'bound state intertwiner' p]| , The 
bound state S-matrix, that is the scattering matrix of the bound state (12) with a 
particle 3, is obtained by the "bootstrap equation" 

.(12) _ ^(12) 



5'(12)3(^(12)3) = 'S'i3(6'i3)S'23(023) 




where we use the usual short notation of matrices acting in the spaces corresponding to 
the particles 1, 2, 3 and (12). 



6 



As an example we consider the sine-Gordon alias massive Thirring model. For u > 1 
(i.e. g < 0) there are no soliton anti-soliton bound states and the S-matrix S{9) is analytic 
in the physical strip < Im^ < tt. For u < 1 there are poles of the S-matrix element 
Sss{0) a.t 9 = Z7r(l — ku) for = 1, ■ ■ ■ < l/i/ which correspond to bound states with 
masses ruk = 2M sin^uk where M is the soliton mass. This mass spectrum coincides 
with the semiclassical breather spectrum if the parameter u is related to the sine- 
Gordon coupling constant (3 as noted above. The pole at 6' = Z7r(l — u) correspond in 
particular to 

soliton + anti-soliton — > lowest breather. 



The "bootstrap equations" yield the breather-soliton S-matrix |62 

SUO) = tie + - z/)) u{9 - - z/)) 

and in a second step the 2-breather S-matrix 162 



sinh 6 + i sin |vr(l + u) 
sinh 9 — i sin |vr(l + u) 



SM = sue + l^7r{l - u)) sue - |^vr(i - u)) = + ' 

smli 9 — t sm iru 



5 Form factors 

Definition 1 For a local operator 0{x) the generalized form factors /jT^ are defined as 

(^1, . . . , ^n) = ( I 0(0) I Pl, . . . , p„ )Z...a^ 

for 9i > ■ ■ ■ > 9n (for other orders of the rapidities they are defined by analytic continua- 
tion). The index ai denotes the type of the particle with momentum Pi. We also use the 
short notations Oa{9) or Oi,,,n{9). 

For the sine Gordon model a = soliton, anti-soliton or breathers. First we restrict our- 
selves to solitonic states. Similar as for the S-matrix 'maximal analyticity' for generalized 
form factors means again that they are meromorphic and all poles in the 'physical strips' 
< Im^^j < TT have a physical interpretation. Together with the usual LSZ-assumptions 
||63| of local quantum field theory the form factor equations can be derived 

(i) Watson's equations: 

0...ij..{- ■ . ,9i,9j, . . .) = 0,„ji„X- ■ ■ 5 Sj, 9i,...) Sij{9i — 9j) 



[ii) Crossing relations: 

^i(Pl|(^(0)|p2,...,Pn)~- = a,a,...a„(^l+^Vr,e2,...,e„) 

= Oa2...a„aii92, ■ ■ ■ ,9n,9i — in) 
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{iii) Recursion relations: 



Res 

9l2=iTT 



12 *-^3...nl 



1^1 ~ S2n ■ ■ ■ 'S'23) 



where C12 is the charge conjugation matrix 
{iv) Bound state form factors equation: 



Res Cl23...n( 

where a is the fusion angle. 
[v) Lorentz invariance: 

Cl...n(^l +U,.. 

where s is the "spin" of O. 



0(12)3...n(^(12),^') 



(12) 
12 



,0„ + m) =e^"Oi...„(ei,...,^„) 



These equations have been proposed by Smirnov |Tl[ as generalizations of equations 
derived in the original articles |T0|, |5^, |^. They have been proven [jl2| by means of LSZ- 
assumptions and 'maximal analyticity'. They hold in this form for bosons; for fermions 
or more generally for anyons there are some additional phase factors. If we write the 
form factors as a formal sum of Feynman graphs and use the additional rule that a line 
changing the 'time' direction also changes a particle to an anti-particle and changes the 
rapidity 9 ^ 9 iin we may depicted the equations (i) — (iv) as in figure |l|. 



(^) 



in 



[III 



[iv] 



— Res 

2l ei2=in 



1 



—= Res 

a/2 Si2=ia 



O 



o 

^conn. 



o 



o 



o 



o 



n 



o 



7^ 




Figure 1: The form factor equations. 



8 



Locality 

It has been proven ||65|0. that the properties (z) — {iv) of the form factors together with 
a general crossing formula imply locality in the form of 



^^{4>\ [0{x),0{y)] |^r = 

for all matrix elements, if a; — y is space like. In the proof one assumes the convergence 
of the sum over all intermediate states. 



Two-particle form factors 

For the two-particle form factors the form factor equations are easily understood. The 
usual assumptions of local quantum field theory yield 

( I 0(0) I Pi, p2)^"/°"* = F {{p, + ±ie) = F {±9^2) 

where the rapidity difference is defined by piP2 = f^i^ cosh ^12. For integrable theories one 
has particle number conservation which implies (for any eigenstate of the two-particle 
S-matrix) 

( I 0(0) I pi, = ( I 0(0) I J92, pir' S (^12) . 

Crossing means 

{p^\0{Q)\p2) = F{nx-e^2) 

where for one-particle states in- and out-states coincide. Therefore Watson's equations 
follow 

F{e) = F{-e)S{e) 

F(iTr-e) = F (in + 9) . 



For general theories Watson's |^ equations only hold below particle production thresh- 
olds. However, for integrable theories there is no particle production and therefore they 
hold for all complex values of 6. It has been shown that these equations together with 
"maximal analyticity" have a unique solution. As an example we write the sine-Gordon 
(alias massive Thirring model) two-soliton form factor [1T0| 



00 







, ^n I dt smh|(l-i/) . „ 9 

Fkn -0) = cosh ^6 exp / , , ^\ , \ sin^ t— . 

^ ' ^ ^ ' t sinh^ cosh^sinht 27r 



^For the case of no bound states this was proven before by Smirnov and Lashkevich m4 
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A formula for generalized form factors 



We are looking for solutions of the form factor equations (z) — {v). For generalized form 
factors for arbitrary numbers of solitons and anti-solitons we make the AnsatzQ 



On 



/^^ dzi--- f^^ dzm h{9, z) p^{9, z) '^gXi., z) 



which transforms the equations (z) — (f ) for the co- vector valued function Oai9) into simple 
equations {i') — {v') for scalar functions Pni^,z) where n is the number of particles. The 
later equations are easily solved (see below). To capture the vectorial structure we use 
the "off-shell" Bethe Ansatz state 2.) (see below). For all integration variables zj 

(j = 1, . . . ,m) the integration contours Cg consists of several pieces (see figure |). The 



o 9r, + in{2i^-l] 

• 9n+ iT^ 

o 9n + ivr(z/ - 1) 



O dn - iTT 



Un — nxu 
9„ — 27rz 



o 9i + nx{2u-l) 

o 92 + i'K{2u-l) 

• 9i+ in 

• 92 + iTT 

o 9i + i-K^u - 1) 
o 92 + iTi{v-l) (T) 9i 



9i — in 



o (^2 - m f - 



*92 — ^TTZ/ 

» 92 - 2m 



*9i — inv 
* 9i - 2TTi 



Figure 2: The integration contour (for the repulsive case v > 1). The bullets belong to 
poles of the integrand resulting from u{9i — uj) (f){9i — Uj) and the small open circles belong 
to poles originating from t{9i — Uj) and r{9i — Uj). 

number of integrations m depends on the charge of the operator q = n — 2m. The scalar 
function h{9_,z) is uniquely determined by the S- matrix 



h{9,z)= H F{9,,)l[l[m-zj) n 

l<i<j<n 1=1 j=l 1<*<J<"1 

with 

. 1 1 

r(z) = — — — a sinh 2; sinh 1/ 



' F{z)F{z + i7r) 

where F{9) is the soliton-soliton form factor above. The dependence of the operator 0{x) 
enters only through the scalar p-functions Pn{^,z) (see below). We consider p-functions 
which satisfy the following conditions: 

''This formula using the so caUed 'ofF-shell Bethe Ansatz' has been introduced by the authors et al. 
. Similar integral representation have been also introduced by Smirnov [jlll . 
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z) is symmetric with respect to the ^'s and the z^s. 



[a') Pni^, z) = {. . . , 9i — 27ri, . . . ,z) and it is a polynomial in e^^^ (j = 1, . . . , m). 



[III 



On + e, On] ~z,z^ = e^) = — p^_2(^, I) + P^'^ ( 

m 
m 



p^(6'i = 6'„ + ivr, ^, On, z, z^ 



where d = {62, ■ ■ ■ , 6'„,-i), 1 = {zi, . . . Zm-i) and x = — (F'(0))^ /vr. The func- 
tions (^) are non- vanishing only for charge less operators and they are indepen- 
dent of the zs. 

{iv') The bound state p-functions are investigated below. 
W) + 1^1 z + n) = e'^'^p^{9,z) where s is the 'spin' of the operator 0{x). 

Again these equations hold in this form for bosons; for fermions or more generally for 
anyons there are some additional phase factors. 

Theorem 1 Let generalized form factors be given by the Ansatz 1^. They satisfy the 
form factor equations (i) — (v) if the p-functions Pni^^z) satisfies the equations {i') — {v'). 



This theorem has bee proven for odd |jT2[ and for even number of solitonic particles. 



The "off-shell Bethe Ansatz" state 

As usual one defines the 'monodromy matrix' 

Tl...n,oiO, 60) = Sio{9i — 9q) S2o{92 — 9o) ■ ■ ■ SnoiOn — ^o) 











1 


2 


n 






acting in the tensor product of the 'quantum space' and the 'auxiliary space' V ® Vq 
with V"^-" = Vl ® • ■ ■ ® K (for the sine-Gordon model which correspond to the quantum 
group slq{2) all Vi = C^). Further one defines as usual the sub-matrices A, B, C, D by 

T (ft A=( Bi,„^ii^z) 

Il...n,0[^,Z)-y c,,„n^e,z) D,U0,Z) 

The "pseudo-vacuum" consists only of solitons 

fli,,,n = s ® ■ ■ ■ ® s . 
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A Bethe Ansatz co-vector in Vi...„ is defined by 



^l...n(^, Z) = ^^l...nCl...„(^, Zi)--- Ci,„ni9, Z.^) ■ 

s s 

s — 



01 



Or 



Zl - 

— s 



The conventional application of the Bethe Ansatz is to solve an eigenvalue problem 
of a spin chain Hamiltonian or of a transfer matrix. The eigenstates are then given by 
a discrete set of Bethe Ansatz vectors where the parameters Zi have to satisfy the Bethe 
Ansatz equations. In the "off-shell Bethe Ansatz" the parameter Zi are summed 
or integrated over as above in our Ansatz for the form factors. For other models which 
correspond to groups or quantum groups of higher rank one has to apply a nested off-shell 



Bethe Ansatz 57 



Examples of p-functions 

Form factors for arbitrary numbers of solitons and anti-solitons are given by the Ansatz 
above which involves the p-functions. For the sine-Gordon alias massive Thiring-model 
there have been proposed the p-functions for several local operators. 

Examples of operators and their p-functions (up to normalizations) 
The fundamental fermi field (charge q = n — 2m = 1) 11121 

(m n 
j=l i=l 

The fundamental breather field (charge q = n — 2m = 0) 



The current = 'ip'y^ijj{x) 



±1 



E^-"'E^"H-E^"E 



The energy momentum tensor (with p, a = +, — ) 
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The cx)-many conserved currents [ES 



(L = ±l,±3,...). 



i=l j=l 

It is easy to check that these p-functions satisfy the equations {i') 
with the quantum numbers of the operators. 



{v') consistently 



Identification of the operators 



Several additional checks |]12|, ^ have been performed to justify the correspondences of 
operators and p-functions. The Feynman graph expansion of the matrix element has been 
compared with the expansion of the exact result given by the integral representation and 
always agreement has been found. For example the 4-particle form factor of the sine 
Gordon <^-field can be calculated in lowest order in g by the Feynman graphs of figure |] 
using Coleman's formula e^^d^'-p = —'^ip'^^il^. The result is 



r 



r 





Pl P2 P3 P4 



Pi P2 P3 Pi 



Figure 3: Feynman graphs 



;o|^(o)|pi, ...,P4) 



m 
ssss 



2nig sinh ^6is sinh ^624, sinh |(6'i2 + 6u] 



+ 0{g' 



(3 ni<j cosh i% 

which agrees with the exact result. Another check is to calculate the eigenvalues of 
charges corresponding to local operators. For example it can be shown |T2[ that with the 
p-functions above the higher charges satisfy the eigenvalue equation 







as is expected from the classical higher conservation laws. 



Breather form factors 

It has been shown it the original article [^] that for particles which do not possess 
backward scattering the generalized form factors can be written as 

a(^i,...,^„)=xf(£) n ^(^^.o (2) 

l<i<j<n 
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where F{6) is the two-particle form factor. It satisfies Watson's equations 



F{e) = F{-9)S{e) = F{27ii 



where S{6) is the diagonal two-particle S-matrix. For the lowest breathers of the sine- 



Gordon model the S-matrix is given above and the two-particle form factor is 



, ' dt cosh it — coshfi + z/)t , 



1-1 



(normalized such that Fi,h{oo) = 1). The K-function K^{9) depends on the operator and 
it satisfies Watson's equations for S = 1. For simple cases the K-functions have been 
proposed in the original article |jlO[- Smirnov used the bound state fusion method 
to obtain a breather form factor formula. For the sinh-Gordon model in several articles 
2T|, E4l H K-functions where proposed for various operators. 



Starting with the general formula (||) for soliton form factors and using the bound state 
fusion method we derive soliton-breather and pure breather form factor formulae. The 
K-function for the case of the lowest breathers turns out to be of the form 



h=0 l„=0 



n 

l<i<j<n 



^ I sm vrz/ 
sinh Qi^ 



(3) 



where the breather p-function /) 
with Q2i-\ = 6i + UttU - v) , 92i 



PsoZ2n(^'^) is obtained from the solitonic one 
U„^^ ,A ^. _ i^TT (1 - (-l)'iz/). In 



this way we [^] obtain the breather p-functions for the fundamental breather field V5(x), 
the energy momentum tensor and the cxo-many conserved currents from the solitonic p- 
functions above. 

However, one may also assume a different point of view. Namely, one can consider 
the form of the K-functions as an Ansatz and look for breather p-functions such that 
the general form factor equations {€) — {v) are satisfied. Doing this one will obtain a 
wider class of p-functions corresponding to operators which are local with respect to the 
breather field, but not necessarily local with respect to solitonic field. For example we 
propose the breather p-function corresponding to the normal ordered exponentials 
of the field : e*"^*^ : [x) for generic real 7 



(/) = AT^^) ]^ with q = exp («y 7^ • 



(4) 



Here and in the following : ■ • • : denotes normal ordering with respect to the physical vac- 
uum which means in particular for the vacuum expectation value ( | : exp i'yip : (x) | ) = 1. 

^For the sinh-Gordon model an analogous representation as (|^) together with this p-function was 
obtained previously |44 by different methods. 
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This breather p-function is not related to a sohtonic p-function of any local operator (at 
least for generic 7). One easily calculates the examples for n = 1,2 explicitly 



[e^^1,(^i,^2) = N^'\q-l/qfFUeu) 

Expanding these relation in powers of 7 one obtains the p-functions of all normal ordered 
powers of the field ip 

Hi) = iVf ) (^B-l)'] with iVl^) = [jYN^f^ (5) 



(since Nn^ is independent on q as we will see below). In particular for = 1 and n = 1, 2 
one calculates 

(0|</.(0)|pi) = 2iV}^) 
(0|<^(0)bi,p2 )^" = 0. 

Note that this breather p-function for ip{x) is not the same as that obtained from the 
solitonic one above using bound state fusion. However, it can be proven [49| that the 
form factors are the same in both cases. To justify the proposal (^ and to calculate the 
normalization constants we investigate the asymptotic behavior of form factors. 



Asymptotic behavior of form factors for :e*^^: 

Let O = : (fi^ : be the normal ordered power of a bosonic field. Set the rapidities as 
9 = \9[, . . . , XO'^, 9'{, . . . , 9"_^ and let A — cx) then the asymptotic behavior of the n- 
boson form factor is 

K=Q ^ ^ 

if the interaction is pure bosonic. This can be proven in any order of perturbation theory 
as follows. The matrix element on the left hand side may be written in terms of Feynman 
graphs as 




9\ + I \9^ ^il ' ' '\9 g I ' ' '\9^ 

where all other graphs not drawn have lines which connect both parts. Weinbergs power 
counting theorem for bosonic Feynman graphs implies that these contributions decrease 
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for A — s> oo as 0{X^e~^) . This behavior is also assumed to hold for the exact form factors 
(the fact is that the 'logarithmic terms' A'^ do not show up for the exact expressions since 
the K-functions are meromorphic in the e^'). Therefore for the exponentials of the boson 
field : e*''"^ : we have the asymptotic behavior 

It is easy to see E9[ that our proposal (H) together with (0) and (Bl) satisfies this 



asymptotic behavior. ^ The asymptotic behavior of other form factors is more complicated 
^ in particular if fermions are involved. 



Normalization of form factors 

The normalization constants are obtained in the various cases by the following observa- 
tions: 

a) The recursion relation {in) relates Nn+2 and Nn- For a typical p-function as that 
for the exponentials of the field this means 

Ni'^=NB2-. Irr-^ in>3). (6) 

b) For a field annihilating a one-particle state the normalization is given by the vacuum 
one-particle matrix element, in particular for the fundamental breather field one has 

{0\!f{0)\p) = 

where Z'^ is the finite wave function renormalization constant. For the sine-Gordon 
field it has been calculated in the original article |iy 



(l + z/)-^exp / —-dt 



sin -^u \ Jo sint 

c) If a local operator is related to an observable like a charge Q = J dxO{x) we use 
the relation 

{p\Q\p) = qip'l p)- 

This may be applied for example to the higher conserved charges. 

d) We use Weinberg's power counting theorem for bosonic Feynman graphs. As dis- 
cussed above this yields in particular the asymptotic behavior for the exponentials 
of the boson field O = : e^^^ : 



Onie.e^, ...) = On-i{e2, ...) + o{e 



-Re 6 



as Re ^1 — i> oo in any order of perturbation theory. This behavior is also assumed 
to hold for the exact form factors. Applying this formula iteratively we obtain from 
relations for the normalization constants. 



^This type of arguments has been also used before 0, p^ . 
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Examples: 

1. For the normalization factors for exponentials of the field given by (^) one uses a) 



and d) ^ |9 



sin TTuFhbd'n' 



2txv 



where the identity Z'^F{iiT)f3'^ sin it u = 8 (ttz/)^ has been used. Note that Nn^is 
independent of q. 

2. For the normalization factors for the field given by (|^) (for = 1) one uses b) which 
imphes 

(0|v?(0) |p) = 2N[^^ = v^. 

This is consistent with our proposals and (^) and justifies the identification 
q = exp (i^7 



6 Some Results 

The quantum sine-Gordon equation 

We start with the local operator : sin •yip : (x) = ^ : (e*''"^ — e"*''"^) : (x). For the exceptional 
value 7 = /9 we find ^ that also : sin Pip : (x) is local. Moreover the quantum 

nip{x) + -:sm(3ip:{x) = (7) 

P 



sine- Gordon field equation 



holds for all matrix elements, if the "bare" mass ^/a is related to the renormalized mass 

byQ 



a = m 



2 TTZ^ 

sin TTU 



where m is the physical mass of the fundamental boson. 

This is a sketch of the proof which uses induction and Liouville's theorem. Con- 
sider the K-functions of the left hand side of (|^) 

Ui) = -E''-j:'-'-Ki:Hi) + (At'(9) - Al""(9)) . 

The results of the previous section imply fi{6) = j^iE) = for Q = ^^^"^ As induction 
assumption we take fn-2{^') = 0. The function is meromorphic in terms of the 

Before such formula was found lH, H by different methods. 
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Xi 



e^' with at most simple poles at Xi = ±Xj since sinh 6ij = (xi + Xj){xi — xj) / (2xiXj) 



The residues of the poles at Xi 



vanish because of the symmetry under the exchange 



of The residues at proportional to /n-2(^") because of the recursion 

relation [iii). Furthermore it can be shown |^ that /n(^) — >■ for Xi — > oo. Therefore 
fn{0 vanishes identically by Liouville's theorem. 

The factor .^'^ in (13) modifies the classical equation and has to be considered as a 
quantum correction. For the sinh-Gordon model an analogous quantum field equation 



has been obtained previously [2^|0. Note that in particular at the 'free fermion point' 
— »• 1 [f3^ — >■ Air) this factor diverges, a phenomenon which is to be expected by short 
distance investigations For fixed bare mass square a and — > 2, 3, 4, . . . the physical 
mass goes to zero. These values of the coupling are known to be specific: 1) the Bethe 
Ansatz vacuum in the language of the massive Thirring model shows phase transitions 
|70| and 2) the model at these points is related |71, 72 , 73] to Baxters RSOS-models which 
correspond to minimal conformal models with central charge c = 1 — 6/(z/(z/ + 1)). 



The trace of the energy momentum tensor 



As a further operator equation we find [^, ^ that the trace of the energy momentum 
tensor satisfies 



-2-^ (l-^)0cos/3y.:(x)-l). 



(9) 



Again this operator equations is to be understood as equations of all its matrix elements. 
The equation is modified compared to the classical one by a quantum correction (1 — 
/S^/Svr). As a consequence of this fact the model will be conformal invariant in the limit 
(5"^ — >■ Stt for fixed bare mass square a. This is related to a Berezinski-Kosterlitz-Thouless 
phase transition. 




2 2 

p = m 




cos Pip 



Figure 4: Bosonic Feynman graphs 

All the results may be checked in perturbation theory by Feynman graph expansions. 
In particular in lowest order the relation between the bare and the renormalized mass (^) 
is given by figure ^ (a). It had already been calculated in the original article |T^. The 
result is 



m 



a 1 - 



1 fl3 



6 V 8 



^It should be obtained from (|7|) by the replacement (3 ig. However the relation between the bare 
and the renormalized mass differs from the analytic continuation of M) by a factor. 
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which agrees with the exact formula above. Similarly we check the quantum corrections 
of the trace of the energy momentum tensor @ by calculating the Feynman graph of 
figure ^ (b) with the result |10 



{p\ :cos/?¥.: (0) - = -P' [^ + ^)+ ^ 

This again agrees with the exact formula above since the usual normalization for the 
energy momentum given by c) implies {p \ T'^|p) = 2m^. 
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